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ABSTRACT 

In  the  present  paper,  an  analytical  theory  will  be  investigated  for  studying  three-layer  slab  waveguide  for  sensor 
applications.  The  sensor  is  a  proposed  nonlinear  optical  waveguide  including  a  linear  thin  film  with  thickness  h  and 
dielectric  constant  Ef  that  is  surrounded  by  a  substrate  with  nonlinear  permittivity  ^and  a  cover  with  nonlinear 
permittivity  e2-  Nonlinear  environments  are  of  Kerr  type.  One  of  the  best  quantities  studied  in  the  research  is  sensor 
sensitivity,  that  is  determined  with  change  in  effective  refractive  index  N  for  changing  in  refractive  index  of  cover 
medium  nc.  Status  of  obtaining  max  sensitivity  for  nonlinear  waveguide  sensors  TE  has  been  investigated  and  compared 
with  linear  sensors.  Guide  layer  has  been  chosen  as  Si3N4and  GaN  and  we  have  applied  silica  as  cover  and  substrate 
environment  with  some  percent  impurity.  Our  proposed  sensor  with  sensitivity  of  Sh=  0.  2798  increases  up  to  105.433% 
and  thickness  of  sensor  has  decrease  of  h=105  nm  up  to  41.666%  compared  to  previous  proposed  sensors  of  researchers. 
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INTRODUCTION 

In  the  recent  two  decades,  optical  waveguide  sensors  have  attracted  much  attention.  In  the  optical  fiber  changes 
process,  various  investigations  have  been  done  that  its  concentrated  has  been  on  suitable  design  of  sensors.  As  offshoot  of 
these  observations,  new  ideas  using  for  the  design  of  sensor  system  that  had  been  led  to  fiber-based  sensor  devices  and 
parts  (Yu,  et  al  2002;  Gholamzadeh,  et  al  2008;  Horvath,  et  al  2003;  Taya,  et  al  2011a;  Abadla,  et  al  2004a). 

Optical  sensors  used  in  modification  of  chemical  measurements  to  optical  properties  such  as  intensity,  phase, 
frequency,....  Chemical  sensors  considered  for  immunity  in  electromagnetic  interference,  mechanical  resistance,  small  size, 
...  and  for  this  reason,  it  can  be  used  in  hazardous  environments  (Niu,  et  al  2011;  El-Agez,  et  al  2011;  Taya,  et  al  2010; 
Taya,  et  al  2012). 

In  the  past  years,  optical  sensors  have  been  used  in  different  cases  for  various  conductions  such  as  biology, 
biochemical  and  biosensors  (Benaissa,  et  al  1998;  Qing,  et  al  1999;  Abadla,  et  al  2004b;  Abadla,  et  al  2003;  Veldhuis, 
et  al  2000).  With  the  use  of  these  sensors  physical  quantities  can  be  measured  such  as  electric  current,  magnetic  field, 
pressure,  heat,  displacement,  pollution  of  sea  waters,  liquid  level,  gamma  and  X  rays. 

Its  further  application  in  medical  can  be  to  detect  the  concentrations  of  certain  chemicals  in  blood,  pharmaceutical 
applications,  cancer  tumors  dosimeter,  identifying  of  deficiencies  within  the  body,  laser  surgery,  dentistry  application, 
liquid  measurement  and  blood. 
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The  essential  part  of  an  optical  waveguide  sensor  consists  of  a  high  refractive  index  wave  guiding  film 
sandwiched  between  a  substrate  and  a  covering.  The  waveguide  sensor  structure  can  be  divided  into  two  configurations. 
The  first  configuration  is  called  a  normal  symmetry  sensor  which  the  substrate  with  a  refractive  index  higher  than  that  of 
the  cover.  The  second  configuration  is  the  substrate  with  a  refractive  index  less  than  that  of  the  cover.  The  waveguide 
sensor  is  called  a  reverse  symmetry  waveguide.  In  general,  the  reverse  symmetric  waveguide  sensors  have  relatively  higher 
sensitivity  (Parriaux,  et  al  2000;  Veldhuis,  et  al  1999;  Skivesen,  et  al  2003;  Horvath,  et  al  2003;  Taya,  et  al  2011b;  Schmitt, 
et  al  2010). 

Effective  refractive  index  (N)  is  one  of  the  most  important  factors  in  the  field  of  optical  sensing  that  is  obtained  by 
equation  N  =  c/vph,  where  c  is  the  speed  of  light  in  vacuum  and  vph  is  the  phase  velocity.  The  change  in  effective 
refractive  index  for  each  cover  evanescent  wave  induced  interactive  changes.  This  effect  is  the  basis  of  optical  sensors. 

When  a  light  beam  with  alpha  angle  was  injected  in  an  end  of  waveguide  in  film  part  and  the  light  is  guided 
within  film  by  total  internal  reflections  in  boundaries  of  film-cover  and  film-substrate  in  around  of  angle  9  and  exited  from 
the  other  end  of  the  waveguide  the  intensity  of  which  was  measured  by  detector  (Abadla,  et  al  2004;  Brioude,  et  al  2000; 
Skivesen,  2005;  Taya,  et  al  2012;  Taya,  et  al  2013;  Taya,  et  al  2011c).  Evanescent  wave  formed  in  each  point  of  integrated 
internal  reflection.  A  small  part  of  evanescent  wave  distributed  into  the  covering  environment  in  comparable  space  with 
guidance  light  wavelength.  This  Evanescent  wave  is  important  in  sensitivity  operating.  When  cover  is  placed  in  the 
evanescent  field,  change  is  observed  in  absorption  or  phase  shift  the  light  propagating  in  waveguide.  This  observed  change 
is  the  sign  of  concentration  or  cover  refractive  index.  This  displacement  between  cover  and  evanescent  field  can  be 
observed  with  spectroscopic  techniques. 

THEORY 

Here,  we  are  primarily  interested  in  examining  the  behavior  of  nonlinear  slab  waveguides  in  sensing  applications. 
In  the  present  study  the  p-polarized  waves  propagate  through  x  axes  are  considered. 

Main  part  of  a  nonlinear  waveguide  sensor  consists  of  a  linear  thin  film  of  thickness  h  and  dielectric  constant  £f. 
This  film  is  sandwiched  between  a  nonlinear  substrate  and  a  nonlinear  covering  layer.  Where  covering  of  nonlinear 
permittivity  (e-l)  and  substrate  of  nonlinear  permittivity  (£2)  are  given  by: 

Ei  =  £c  +  CCc|Eyl|2.  (1) 

£2  =  £s  +  ccs|Ey3|2.  (2) 

Where  ac  and  as  are  the  nonlinearity  coefficients  of  covering  and  substrate,  respectively,  £c  and  £s  are  the  linear 
parts  of  the  permittivities  and  Eyl  and  Ey3  are  the  TE  fields  in  covering  and  substrate,  respectively.  The  electric  field  in 
each  layer  has  the  following  solution  for  TE  modes  for  occ  >  0  and  ocs  >  0: 

Eyi  =  J—  qcsech(k0qc(z  -  zc)),     z  >  h  (3) 

Ey2  =  Acos(k0qfz)  +  Bsin(k0qfz),     0  <  z  <  h  (4) 
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Eys  =    f"  qssech(k0qs(zs  -  z)),      z  <  0.  (5) 

Where  qc  =  -^N2  —  £c  ,qf  =  -Jsf  —  N2,  qs  =  -^N2  —  £s,  N  is  the  effective  refractive  index,  constants 
A  and  B  represent  the  amplitude  of  the  waves,  and  zc  and  zs  are  constants  related  to  the  field  distribution  in  the  covering 
medium  and  the  substrate,  respectively. 

Continuity  of  Ey  and  — —  gives  the  dispersion  relation 

tan(k0qfh)  =  Q^f  tanh  Cc+qsqf  tanh  Cs^  (  } 

v   um   J        q2-qcqs  tanh  Cs  tanh  Cc  v  J 

Where  k0  is  the  free  space  wave  number,  Cc  =  k0qc(h  —  zc),  and  Cs  =  k0qszs  which  are  called  the  covering  film 
interface  nonlinearity  and  the  substrate  film  interface  nonlinearity  respectively. 

Where  ac  and  as  are  two  asymmetry  parameters  and  Xc  and  Xs  are  two  normalized  variables,  where  ac,  as,  Xc,  Xs 
are  given  by: 

a*.  —     .  ar  —     .  Xc  —     Xr  —    .  (7) 

£f         Ef  qr  qr 

Xc  and  Xs  are  linked  by: 

2  =  (l-ac)(l,Xi)_1 
(l-as) 

The  effective  refractive  index  can  be  written  in  terms  of  a9  and  X,.  as: 


N  =  V^J^f-  (9) 
In  terms  of  Xc  and  Xs  Eq.  (6)  can  be  written  as 

k0qfh  —  arcta  n(Xc  tanh  Cc)  —  arctan  (Xs  tanh  Cs)  —  rrm  =  0.  (10) 
Where  m  =  0,1, ...  is  the  mode  order. 

From  the  dispersion  relation  given  by  Eq.  (10),  we  derive  the  sensor  sensitivity  Sh,  i.e.,  the  basic  sensing  principle 
of  the  planar  dielectric  waveguide  sensor  is  to  measure  changes  in  N  due  to  changes  in  nc.  Differentiating  Eq.  (10)  with 

(dnr 


respect  to  N  and  calculating  Sh  as  C  n<:/gisj)    we  obtain: 


/a7Jl+X|  (Hc+tanhCc) 

(11) 


XcJac+X£  (1+X^tanh2  CC)(ATE  +  GS+GC) 

Where 


Hc  =  k0(h  -  zc)XcA/^  If^l  (1  -  tanh2  Cc).  (12) 

_  Hc+tanhCc  (1+Xg)  -  . 

c       Xc(l+X^  tanh2  Cc)  '  ^  ' 
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Hs+tanhCs(l+X^ 

s       Xs(l+x|tanh2Cs)  v  J 


Hs  =  k0(zs)XsV^         (1  -  tanh2  Cs).  (15) 


ATE  =  arctan  (Xs  tanh  Cs)  +  arctan  (Xc  tanh  Cc)  +  nm.  (16) 

The  sensitivity  of  the  proposed  nonlinear  sensor  is  given  by  Eq.  (11)  which  strongly  depends  on  the  waveguide 
structure. 

One  of  the  most  important  quantities  is  the  power  flow  in  different  layers  of  the  slab  waveguide.  The  sensitivity  of 
the  sensor  is  critically  dependent  on  the  fraction  of  total  power  propagating  in  the  covering  medium.  The  energy  flux 
per  unit  length  is  given  by 

P  =  J_+jExHdz  =  Ps+Pf+Pc.  (17) 
Where  Ps  is  the  power  flow  in  nonlinear  substrate  and  Pf  is  the  power  flow  in  linear  film. 
The  fraction  of  total  power  flowing  in  the  nonlinear  covering  is: 

_?c_  _   (18) 


Ptotal       Xc-     X^sech^Cs     Xs  ' 
acJ1         2as  as'2 

Where 

r  =  k0qfhx+  +  isin(2k„qfh)x_  +  Xs  tanh  Csj3 .  (19) 

ji  =  1  -  tanh  Cc.  (20) 

j2  =  l-tanhCs.  (21) 

j3  =  1  -  cos(2k„qfh).  (22) 

x_  =  1  -Xs2tanh2cs.  (23) 

x+  =  1  +Xs2tanh2cs.  (24) 

REPRESENTATION  AND  DISCUSSIONS 

In  our  calculations,  we  will  assume  the  guiding  layer  to  be  Si3N4  (nf  =  2)  and  GaN  (nf  =  2.34949),  the  free  space 
wavelength  will  be  X  —  1550  nm,  tanhCc  =  0.6  and  tanhCs  =  0.7.  Only  the  fundamental  mode  (m  =  0)  will  be 
considered  since  it  has  the  highest  sensitivity. 

Substratum  and  covering  medium  are  considered  silica  with  different  impurity  percent.  It  is  necessary  to  mention 
that  silica's  linear  and  nonlinear  dielectric  constant  is  respectively  £  =  2.13  and  a  —  7.04  X  10~20. 

In  figure  (1),  proposed  nonlinear  sensor  sensitivity  is  drawn  based  on  guidance  layer  thickness  h  that  conventional 
symmetric  waveguide  (i.e  ns  >  nc)  is  considered.  Covering  medium  for  this  diagram  is  silica  with  £c  =  2.13  and 
substrate  environment  has  been  regarded  as  silica  with  three  different  kinds  of  impurity  which  are  0.09%,  0.9%  and 
10%  dielectric  constant  impurity  consist  of  £s  =  2.311917,2.14917  and  2.343  and  diagrams  of  each  one  are  drawn 

|  Index  Copernicus  Value:  3.0  -  Articles  can  be  sent  to  editor@impactjournaIs.us  | 


Theoretical  Analysis  of  Nonlinear  Optical  Waveguide  Sensors:  TE  Case 


129 


separately  and  then  they  are  compared.  As  it  is  seen  in  the  figure,  substrate  with  £s=  2.311917  has  the  best  sensitivity  and 
minimum  thickness  of  guidance  layer  h.  It  should  be  noted  that  nf  =  2  has  been  considered  and  then  diagram  for 
nf  =2.34949  is  drawn  in  figure  (2)  and  both  figures  are  compared. 


2.343 


0  >C4»»3fr)444«4a4«4~WM* 

h(run) 


Figure  2:  Sensitivity  with  the  Film  Thickness  for 
nf  =  2.  34949,  £c  =  2. 13,  es  =  2. 14917,  es  =  2. 131917,  es  =  2.  343 

The  sensitivity  of  the  proposed  sensor  was  plotted  with  the  film  thickness  for  the  different  asymmetry  parameters 
£s  in  Figure  1  and  2.  The  figures  show  that  the  sensitivity  decreases  with  increasing  £s.  This  is  a  normal  behavior. 
The  sensitivity  operation  is  depending  on  the  evanescent  optical  field  extending  from  the  thin  guiding  film  into  the 
covering  medium.  To  obtain  high  sensitivity,  it  is  essential  to  get  an  amount  of  spreaded  optical  power  in  the  cover. 
Therefore,  increasing  £c  and  decreasing  £s  is  dependent  on  flowed  power  fraction  in  covering. 

Comparing  figures  (1)  and  (2),  it  is  clear  that  in  nf  =2.34949,  sensitivity  has  been  increased  and  thickness  of 
guidance  layer  has  been  decreased.  We  would  like  to  compare  linear  and  nonlinear  sensors  that  have  been  shown  in 
figure  (3).  In  linear  state,  tanhCs  =  1.  We  can  see  that  sensitivity  in  nonlinear  state  is  better  than  the  linear  one  and 
thickness  of  sensor  guidance  layer  is  thinner  in  nonlinear. 


h(nm) 

Figure  3:  Sensitivity  with  the  Film  Thickness  for  nf  =  2,ec  =  2.13,es  =  2.  343 
for  the  Proposed  Nonlinear  Sensor  (Solid  Line)  and  a  Linear  Sensor  (Dotted  Line) 
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One  of  the  important  works  done  in  this  area  is  the  comparison  between  suggested  sensor  diagrams  and  previous 
research  works  that  have  been  shown  in  figure  (4). 


Figure  4:  Sensitivity  with  the  Film  Thickness  for  nf  =  2,  ec  =  2. 13,  es  =  2/343  for  the  Proposed  Nonlinear 
Sensor  (Solid  Line)  and  nf  =  2,  ec  =  1. 8,  es  =  2. 6  for  the  Previous  Research  Works  (Dotted  Line) 

In  figure  (5),  there  is  a  comparison  between  reverse  symmetric  sensor  (ns  <  nc)  and  normal  symmetric  sensor. 
Maximum  value  of  sensitivity  near  the  cut  thickness  decreases  with  increase  in  guidance  layer  thickness  of  sensor. 
This  situation  is  attributed  to  power  considerations.  When  guidance  layer  thickness  is  near  cut  thickness,  effective 
refraction  index  N  is  near  covering  coefficient  nc,  the  penetration  depth  of  the  evanescent  field  into  the  cover  medium 
becomes  infinite  and  total  power  frequently  flows  in  covering  such  that  guidance  layer  thickness  of  reverse  symmetric 
waveguide  increases. 


One  of  the  cases  that  is  considered,  is  sensitivity  diagram  based  on  tanh  Cc.  As  shown  in  figure  (6),  when  tanh  Cc 
goes  toward  1,  the  sensitivity  decreases.  This  is  the  prospected  subject  because  number  1  is  marker  of  linear  sensors. 
According  to  previous  diagrams,  nonlinear  sensors  have  higher  sensitivity  comparing  to  linear  sensors. 
This  diagram  supports  previous  diagrams. 


h(iun) 


hCnm) 


Figure  5:  Sensitivity  with  the  Film  Thickness  for 
=  2. 131917,  ec  =  2. 130213,  ec  =  2. 14917,  ec  =  2.  343 


tanh  Cc 


Figure  6:  Sensitivity  with  the  Film  Thickness  for 
es  =  2.  343,  ec  =  2. 2791,  ec  =  2.  3856,  ec  =  2. 4495,  h  =  0. 1  um 
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Eys  =    f"  qssech(k0qs(zs  -  z)),      z  <  0.  (5) 

Where  qc  =  -^N2  —  £c  ,qf  =  -Jsf  —  N2,  qs  =  -^N2  —  £s,  N  is  the  effective  refractive  index,  constants 
A  and  B  represent  the  amplitude  of  the  waves,  and  zc  and  zs  are  constants  related  to  the  field  distribution  in  the  covering 
medium  and  the  substrate,  respectively. 

Continuity  of  Ey  and  — —  gives  the  dispersion  relation 

tan(k0qfh)  =  Q^f  tanh  Cc+qsqf  tanh  Cs^  (  } 

v   um   J        q2-qcqs  tanh  Cs  tanh  Cc  v  J 

Where  k0  is  the  free  space  wave  number,  Cc  =  k0qc(h  —  zc),  and  Cs  =  k0qszs  which  are  called  the  covering  film 
interface  nonlinearity  and  the  substrate  film  interface  nonlinearity  respectively. 

Where  ac  and  as  are  two  asymmetry  parameters  and  Xc  and  Xs  are  two  normalized  variables,  where  ac,  as,  Xc,  Xs 
are  given  by: 

a*.  —     .  ar  —     .  Xc  —     Xr  —    .  (7) 

£f         Ef  qr  qr 

Xc  and  Xs  are  linked  by: 

2  =  (l-ac)(l,Xi)_1 
(l-as) 

The  effective  refractive  index  can  be  written  in  terms  of  a9  and  X,.  as: 


N  =  V^J^f-  (9) 
In  terms  of  Xc  and  Xs  Eq.  (6)  can  be  written  as 

k0qfh  —  arcta  n(Xc  tanh  Cc)  —  arctan  (Xs  tanh  Cs)  —  rrm  =  0.  (10) 
Where  m  =  0,1, ...  is  the  mode  order. 

From  the  dispersion  relation  given  by  Eq.  (10),  we  derive  the  sensor  sensitivity  Sh,  i.e.,  the  basic  sensing  principle 
of  the  planar  dielectric  waveguide  sensor  is  to  measure  changes  in  N  due  to  changes  in  nc.  Differentiating  Eq.  (10)  with 

(dnr 


respect  to  N  and  calculating  Sh  as  C  n<:/gisj)    we  obtain: 


/a7Jl+X|  (Hc+tanhCc) 

(11) 


XcJac+X£  (1+X^tanh2  CC)(ATE  +  GS+GC) 

Where 


Hc  =  k0(h  -  zc)XcA/^  If^l  (1  -  tanh2  Cc).  (12) 

_  Hc+tanhCc  (1+Xg)  -  . 

c       Xc(l+X^  tanh2  Cc)  '  ^  ' 
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(dnr 


respect  to  N  and  calculating  Sh  as  C  n<:/gisj)    we  obtain: 


/a7Jl+X|  (Hc+tanhCc) 

(11) 


XcJac+X£  (1+X^tanh2  CC)(ATE  +  GS+GC) 

Where 


Hc  =  k0(h  -  zc)XcA/^  If^l  (1  -  tanh2  Cc).  (12) 

_  Hc+tanhCc  (1+Xg)  -  . 

c       Xc(l+X^  tanh2  Cc)  '  ^  ' 
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